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Abstract 


Neutron radiographic imaging is combined with locally resolved current density measurements to study the effects of local water content on the 
performance of the corresponding electrochemical active area in an operating PEM fuel cell. Liquid water agglomerates are detected, quantified 
and correlated with the activity of the respective area. At low currents, depletion of the reactant gas leads to a decreasing performance along the 
anodic flowfield channel. At high currents, an optimum humidification is reached in the central part of the fuel cell; close to the inlets respectively 
outlets, flooding and drying can be observed concurrently and cause a non-uniform current density distribution across the reactive area. The fast 
response of the local performance on water droplets migrating in the gas channel is tracked by short-term imaging taking place on a timescale of 


several seconds. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


The worldwide energy demand and the rapid approach of 
“peak oil’, whichis expected to be reached within this decade and 
results in an overall decreasing production of crude oil, require 
a drastic change towards highly efficient energy conversion sys- 
tems and regenerative sources. Fuel cells bear the potential for 
applications as automotive as well as residential power supplies. 
Among several types of fuel cells, low temperature polymer 
electrolyte membrane fuel cells (PEMFC) are considered for a 
wide range from mobile to portable as well as stationary appli- 
cations. Besides financial aspects such as the overall price the 
long-term stability is still an issue, which has to be solved prop- 
erly to achieve a broad market penetration of these devices. A 
key task to clarify the aforementioned technical problems is a 
proper water management of the system; a minimum humidity 
inside the cell is necessary for the membrane as a dry polymer 


* Corresponding author. Tel.: +49 731 9530 403; fax: +49 731 9530 666. 
E-mail address: christoph.hartnig @ zsw-bw.de (Ch. Hartnig). 


0378-7753/$ — see front matter © 2007 Elsevier B.V. All rights reserved. 
doi: 10.1016/j.jpowsour.2007.08.058 


undergoes structural changes resulting in a vanishing protonic 
conductivity. On the other hand, excess liquid water in the gas 
diffusion layer (GDL) or the flow field channels hinders the gas 
transport and therefore the supply of reactants to the catalyst 
layer. 

During the last years, different experimental techniques were 
applied to study the transport of liquid water inside the flow 
field: several groups employed transparent components to mon- 
itor optically the formation of liquid water. Weng et al. [1,2] and 
Ma and coworkers [3] observed flooding effects on the cathode 
of PEM fuel cells and deduced from these observations impacts 
on the overall performance of the cell. This method has also been 
applied to direct methanol fuel cells by Scott and coworkers [4] 
and Zhao and coworkers [5] targeting for an optimized carbon 
dioxide management, which is one of the key factors for a highly 
efficient DMFC. 

Further experiments extended the optical detection of liq- 
uid water by spatially resolved current density measurements. 
Hebling and coworkers [6,7] evaluated simultaneously current, 
temperature and water distribution inside the cell by means of 
transparent end plates. A further improvement was reached by an 
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additional measurement of the local electrical impedance spec- 
tra [8]. Spatially resolved information of the current density has 
been reported by means of different techniques, either a seg- 
mented flow field [9], respectively collector plate [10], or via an 
array of shunt resistors [11] embedded in the current collector 
plate. By either approach the water distribution and transport 
inside the PEFC can be studied in detail. 

Common to the aforementioned techniques is the strong dis- 
turbance of the fuel cell setup with regard to electrical or thermal 
conductivity. The first attempts to study an unmodified and there- 
fore undisturbed system were achieved by means of neutron 
radiography. This method is based on the high attenuation coef- 
ficient of hydrogen compared to one of the most metals; the 
neutron beam can penetrate the metallic housing of the cell 
as well as the graphitic or metallic flowfields almost unatten- 
uated whereas liquid water leads to a strong absorption of the 
beam which allows for a detection of liquid water down to sev- 
eral nanoliter with a spatial resolution of down to 100 wm. This 
method was engaged by a number of groups [12-15] to detect 
liquid water in an unmodified fuel cell. First studies of material 
variations proved the method as suitable tool for target-oriented 
material development [16]. 

Recently, Scherer and coworkers combined the direct 
observation of liquid water with locally resolved impedance 
measurements and demonstrated that in a fuel cell driven in 
co-flow mode drying and flooding might be observed simul- 
taneously, depending on the respective section of the PEMFC 
[17]. 

In this article we report combined neutron radiography stud- 
ies with locally resolved current density measurements; this 
combination allows for a direct observation of liquid water and 
includes information on the local performance of the cell at the 
time of imaging. The technical details and the neutron radio- 
graphic setup are explained in the following section; a detailed 
analysis and the correlation between the amount of water and 
the local performance are explained in Section 3. In the conclu- 
sion, attempts towards a more uniform current distribution will 
be sketched. 


2. Instrumentation and experiment 
2.1. Radiography setup 


The radiography experiments were performed at the neutron 
tomography instrument CONRAD/V7 at Hahn-Meitner Institute 
(HMI). It is located in the neutron guide hall of the research reac- 
tor BER-2 where low-energy neutrons (so-called cold neutrons) 
are used for experimental purposes. The radiography instrument 
is based on a pinhole geometry with a small variable aperture 
placed at the end of the neutron guide followed by a converged 
collimator with a length of 5m toward the sample position. In 
this way, a well-collimated neutron beam is directed to the sam- 
ple. The transmitted beam is detected behind the sample by a 
position sensitive detector. The main part of the detector sys- 
tem is the 16-bit low-noise CCD camera (Andor DW436N with 
2048 x 2048 pixel), which is focused by a lens system on a 
neutron sensitive scintillator screen. To prevent any radiation 


damage of the camera, the image is deflected on a side under 
90° by a mirror. The sample environment is equipped with a 
translation table, which allows to move the sample to the side 
and to take images of the unperturbed beam. This image is used 
further for normalization of the images with a sample in the 
beam. The normalization procedure helps to transform the gray 
levels from the radiography image to transmission values, which 
are the base for further quantitative analysis. 


2.2. Fuel cell setup 


A single cell setup has been used for the measurements. On 
the anodic and cathodic electrodes, a threefold serpentine flow 
field with 1 mm wide channels and ribs and an active area of 
10cm x 10cm was applied. The flowfields were machined into 
graphite composite plates (SGL Carbon) with cooling flowfields 
on either electrode to ensure a proper thermostatization. GORE 
PRIMEA 5620 membrane electrode assemblies and SGL Car- 
bon 10 BB material were employed in the fuel cell setup. The 
locally resolved current density was measured by means of shunt 
resistors embedded in a segmented multi-layer printed circuit 
board with a segment size of 4cm*; a more detailed descrip- 
tion is given by Lehnert and coworkers [18]. The circuit board 
is attached in between the flow field and the end plate so that 
the local influence on the thermal and electrical conductivity 
of the components and the water development and transport in 
the cell is completely eliminated. Thus, this setup turns further 
modification or interventions of the cell redundant. 

The fuel cell was fed in counter flow mode with hydro- 
gen and humidified air (dew point 25°C) with an anodic and 
cathodic utilization of 80%, respectively 25% at a cell tempera- 
ture of 60°C; the flowfield layout as well as the flow directions 
of the reactant gases together with the matrix for the current 
density measurement are given in Fig. 1. The overall current 
density ig of the galvanostatic measurements was varied in steps 
of 200 mA cm~? from ip = 700 mA cm~? to ig = 100 mA cm~?. 
Here, special attention was payed to the two outmost values 
(700, respectively 100 mA cm7°?), as at high currents a high 
amount of water is produced whereas at low currents and the 
applied stoichiometry the removal of excess water is limited. 
Both conditions might lead to accumulation of excess water in 
the flowfields. At a current density of ip =700 mA cm, the 
potential could not be kept constant which gives a first hint on 
such effects. 

Due to the epoxy substrate of the circuit board, a highly inten- 
sive neutron beam was necessary to transmit the plate and to 
achieve a good contrast for the small amount of water present in 
the GDL. Therefore, the maximum aperture of 3 cm was setup 
at the neutron radiography facility, which gives a spatial reso- 
lution of 500 um at the sample position and a neutron flux of 
2 x 10’ neutrons cm7? s. An exposure time of 0.2s per image 
with a read out time of 1.8 s was applied to perform one radiog- 
raphy for every measured current density distribution. This way 
the time interval between two images in the image sequence 
taken in the experiment was 2.0 s in total. The timers of the cir- 
cuit board and the radiography instrument were exactly aligned 
allowing to correlate the visual information of the water distri- 
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Fig. 1. Threefold serpentine flowfield layout and directions of reactant gas flow; left: anode, right: cathode. The segments of the printed circuit board are displayed; 
colors and numbering refer to the explanations given in the text (For interpretation of the references to color in this figure legend, the reader is referred to the web 


version of the article.). 


bution in the cell with the two-dimensional distribution of the 
power output of the cell. 


3. Results 


The data processing was performed by a standard proce- 
dure: At first, the dark current signal was subtracted from the 
images, which were normalized further by using the unper- 
turbed beam profile. In order to get quantitative information 
on changes of the water distribution in the cell a new renor- 
malization of the image set was performed. Quotient images 
were calculated using the initial dry cell as a reference image 
for further division of the images taken at different operational 
conditions. In this way, only the contribution of the water distri- 
bution can be seen whereas all other details like cell geometry 
and beam inhomogeneity are removed in the resulting images. In 
Fig. 2a, an image as obtained from the experiment is displayed. 
The normalization of this image results in Fig. 2b where only 


the water is visible and all other components from the fuel cell 
are eliminated. In order to illustrate the current density distri- 
bution over the active area the respective segment fields of the 
circuit board have been colored according to the current density 
level and were added as overlays to the normalized radiographic 
images. 

Long-term effects affecting the liquid water distribution were 
separated from quick responses of the cell performance by tak- 
ing the average of ten radiographic images as well as the average 
of the current density. Thereby, contributions taking place on a 
short time scale (covered by single images) and effects of small 
droplets inside the gas diffusion layer were neglected; these 
effects will be discussed further below in more detail. Aver- 
aging the images results in an integral presentation of changes 
taking place within 20 s. 

In Fig. 3a two normalized radiographic images with their 
respective current density distributions at an overall current 
density of iọ = 100 mA cm~? are displayed. The liquid water is 


(b) 


Fig. 2. (a) Radiograph as obtained from the experiment; horizontal lines denote the flowfield channels, vertical lines the contacts of the circuit board and the cooling 
flowfield. (b) Normalized image: liquid water can be identified as black lines and spots in the flowfield channels and the gas diffusion layer. 
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Fig. 3. (a) Water distribution and corresponding current density at cell current /= 100 mA cm7 
(arbitrary units) and respective current density of the segments as denoted in Fig. 1. 


evenly distributed over the whole active area with some (par- 
tially) flooded channels, which can be identified by dark lines. 
As determined from the flow directions of the reactant gases, 
the flooding occurs predominantly on the anode side, which was 
observed already by other groups [12]. At first glance, the cur- 
rent density distribution is, similar to the water distribution, quite 
uniform with its maximum density in the center of the cell; this 
maximum is caused by an interplay of several factors such as 
flow dynamics, cell design and electric contact. In the center, 
almost laminar flow conditions can be reached leading to a per- 
fect supply of the reactive area. Due to the design of the end 
plates, the compression and therefore the electric conductivity 
reaches its optimum in the center of the cell. Reduced current 
densities in the lower part result from a depletion of the reactant 
gas hydrogen, which is fed on the ‘top’ of the flowfield with 
the outlet at the bottom. As the anodic utilization ratio is set 
to ua = 80%, the remaining concentration close to the outlet is 
not sufficient to reach comparable current densities as in earlier 
stages of the flowfield. 

The total volume of the liquid water agglomerates and its 
correlation with the current density is tracked in more detail in 
Fig. 3b. Relative values along one row are considered to track 
changes along the reactant pathway. The hereby used colors cor- 
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respond to the respective segment rows in Fig. 1a; the numbering 
of the plots denotes the segment column in Fig. 1a. Exemplified 
two columns are chosen which represent the whole system ade- 
quately. The variations in the local current densities are at this 
current load small between the segments along the vertical axis. 
Differences lie within a range of 0.05 A cm, across the active 
area the densities are in between 0.13 and 0.2 Acm7?. As evident 
from the neutron radiographies the uniform water distribution in 
the border segments of the active area is reflected in the quanti- 
fied water volumes which all lie within a 10% range (column 5). 
In the central column, the water distribution is not as uniform 
any more with a clear trend towards increased amounts of liq- 
uid water in the lower fields. At the given operating conditions, 
however, a non-uniform water distribution does not affect the 
performance of the cell and the measured current densities are 
within the same range. It is remarkable that the lowest amount of 
liquid water does not correspond to the highest current density, 
which one might expect from a humidification point of view and 
gives a hint on a possible drying of the membrane. 

At an average current density of iọ=500 mA cm~? the 
water distribution is not as uniform as at lower currents; fur- 
thermore, the maximum of the current density distribution is 
shifted slightly upwards compared to operating conditions with 
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io = 100 mA cm7? (see Fig. 4a). In the lower part of the flow- 
field, an increased amount of liquid water can be observed; in 
the lowest section of the fuel cell the liquid water content is 
decreased again as at the cathodic inlet the back-diffusion is 
reduced due to the minor humidification of the reactant gas. The 
quantification of the water agglomerates in Fig. 4b corroborates 
this trend, where the same segments as in Fig. 3b are chosen. 
From a current density of ig = 500 mA cm~? onwards, the impor- 
tance of a well-balanced water management becomes evident as 
can be easily identified in Fig. 4b. In segment column 5 (as 
denoted in Fig. 1b) the respective current density distribution 
is not as close as it was the case at lower currents (compare to 
Fig. 3b). The water distribution also shows a clear trend that is 
related to the current output: neither the lowest nor the highest 
water content lead to the best performance but an intermediate 
humidity gives the highest current output. From a certain water 
content onwards the power outage decreases in the lower part 
of the cell and the combination of excessive liquid water and 
depletion of the reactant gas leads to a reduced current density 
compared to the central segments. In the topmost segment, the 
amount of water is apparently not sufficient to fully humidify 
the membrane and the performance is also reduced compared 
to the segments below. In the central segments (column 3), the 
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effects of the misbalanced water management is even intensi- 
fied (Fig. 4b); the lowest and the highest water amounts are both 
correlated to reduced current densities. 

From ig =500 to io =700 mA cm7? the current density dis- 
tribution becomes even more diverse (Fig. 5a). The flowfield 
channels in the lower third of the active area are more or less 
flooded leading to a decreased performance in the area under the 
respective segments. The time dependency of the current den- 
sities and of the water content are not as uniform as at lower 
currents. In the center, the distribution is even more irregular as 
in the border area; the maximum and minimum of the current 
density differ by almost 100% (Fig. 5b). The fluctuating water 
content in those segments is related to a noisy current density 
distribution as function of time. Larger water drops ‘rushing’ 
through the channels cause irregularities of the supply and lead 
in turn to changes in the current output. A detailed examination 
of the short-term effects of water droplets migrating through 
the flowfield channels will be performed in the following sec- 
tion. In the topmost segments, where effects of the insufficient 
humidification are easily identified by means of the water quan- 
tification curve, the time dependency of the current density is 
quite regular as there are no perturbations stemming from water 
droplets. 
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Fig. 4. (a) Water distribution and corresponding current density at cell current density ip =500mA cm~. The time distance between the two images is 820s. (b) 
Water content (arbitrary units) and respective current density of the segments as denoted in Fig. 1. 


Ch. Hartnig et al. / Journal of Power Sources 176 (2008) 452—459 


(a) 


0.4 


4 4 
18600 18800 


time (s) 


4 
18400 


current density (Acm-?) water volume (a.u.) 


457 


12000 


enereetne 


0-0-0-0-0-0 9-0 9-9, 0-0-0 A 


8000 


4000 


Se oS 
> o wo 


18400 18600 


time (s) 


18800 


current density (Acm?) water volume (a.u.) 


Fig. 5. (a) Water distribution and corresponding current density at the overall current density ip = 700 mA cm~?. The time distance between the two images is 800s. 
(b) Water content (arbitrary units) and respective current density of the segments as denoted in Fig. 1. 


The short-term response of the local performance on water 
droplets migrating through the flowfield channels is studied in 
more detail on basis of single images instead of averaged water 
contents, which are used to describe the long-term behavior. In 
Figs. 6 and 7, two examples of fast activity responses in a distinct 
area on a water droplet passing through an adjacent gas channel 
are given. As demonstrated before, at an average current density 
of ig = 700 mA cm~°, the local current density is not as uniform 
as at lower power requests. Not only in the lower third but also 
in the horizontal midfield of the considered segments several 
channels are at least partially flooded causing an undersupply of 
the reactive area. In the left image of Fig. 6a, the overall water 
content in the highlighted segment is high due to up to five 
flooded channels resulting in a reduced current density; after 
10s (right image of Fig. 6a), the water is transported away from 
this position (and away of this segment) and the normal supply 
of reactant gases is restored. The local current density returns 
back to its original (higher) value as is also illustrated by the 
color change of this segment. The reduction of liquid water in 
this field is additionally monitored by means of the quantified 
water content. Coinciding with the decreasing water content at 
t7 18,350s, the current density is increasing abruptly. These 
sudden changes in the water content and the current density 
(spikes in Fig. 6b) typically result from small pressure changes 


caused by blocked channels, which are emptied if the pressure 
drop becomes too large. 

In Fig. 7, the opposite situation compared to Fig. 6 is 
described, the immediate response of the local performance on 
a droplet ‘entering’ the segment. With the droplet becoming vis- 
ible in the time dependent water content, the current density of 
this segment is decreasing as well. Thereby, the incoming droplet 
blocks the channel and the supply of reactant gases is stopped. 
Effects like cross-flow from a neighboring channel as well as 
cross-diffusion within the gas diffusion layer avoid a complete 
blockage of this area, however, the current density decreases 
by 10-20%. The droplet is not immediately transported away 
but the water content decreases slowly, indicating either a slow 
evaporation of the liquid water or a slow transport of water inside 
the gas diffusion layer, which accumulates inside if there is no 
other way to enter the blocked channel. As images were taken 
every 2s, a time uncertainty is inherent in these images, which 
might camouflage a slight delay of the performance response on 
the blocked channel. 


4. Discussion 


Combined neutron radiographic imaging and locally resolved 
current density measurements were used to investigate the water 
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Fig. 6. Short-term response of the local current on changes of the water content in the channel. (a) Neutron radiographic images taken within an interval of 10s to 
render migrating droplets in the channel. (b) Water content and current of the respective segment; after the disappearance of the droplet the current density increases 
again. 
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Fig. 7. Short-term response of the local current on changes of the water content in the channel. Left: considered segment of the electrochemically active area. Right: 
water content and current density in the respective segment at an overall current density of ig = 500 mA cm~?; the appearance of the droplet leads to a sharp decrease 
of the local current density which recovers slowly when the overall water content decreases again. 
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management in an operating fuel cell equipped with a cur- 
rent density measurement board. The reactant gases were fed 
in counter flow mode with humidification on the cathode side. 
This setup should allow for a uniform distribution of water with 
the anode inlet next to the cathode outlet so that effects like 
back-diffusion can support the equal distribution of liquid water. 

At low power requirements (average current densities in the 
range of 100 mA cm7?) the water as well as the local current 
density distribution is quite uniform over the electrochemically 
active area. Some of the flowfield channels are filled with water, 
whereas an undersupply of the reactive area is not very likely 
since the distribution of current densities is very narrow, bearing 
only minor deviations. In contrary, most probably these agglom- 
erates do not block the entire channel cross-section but fill the 
channels partially leaving enough space for reactant gases. With 
increasing performance requirements (io = 500 mA cm~7) two 
effects can be directly observed: the maximum of the local cur- 
rent density shifts closer towards the anode inlet and in the region 
close to the anode outlet performance losses were observed. 
The latter ones are caused by a combination of reactant gas 
depletion and first stage flooding effects. Most of the liquid 
water is observed in this area, hindering the reactant supply 
and leading to reduced local current densities. The minimum 
water content at this current is observed in the upper part of 
the cell; the performance of these segments is somewhat lower 
than in the central parts, indicating an insufficient humidifica- 
tion of the membrane. This effect becomes more prominent 
at higher currents (i9 = 700 mA cm~?) where a further perfor- 
mance decrease in the upper area is observed. Similar results 
were reported recently by Scherer and coworkers using locally 
resolved impedance spectroscopy [17]. Excess water on the 
anode agglomerates in the lower part of the flowfield where 
single channels of the flowfield are completely blocked and the 
performance in this part is drastically reduced. A more balanced 
water management might be achieved in such cases with thinner 
membranes, enhancing effects like back diffusion to equalize 
excess water on the one side and drying on the cathode side 
[19]. 

Besides long-term changes of the water distribution and the 
correlation with the overall performance, the short-term influ- 
ence of single water droplets in the flowfield channels on the 
local current density has been studied. Small agglomerates are 
able to block at least partially the channel cross-section, caus- 
ing a decrease of the performance when two adjacent channels 
are affected. As soon as these droplets are removed, the perfor- 
mance recovers and the current density reaches again the liquid 
free value. Effects of water droplets migrating through the chan- 
nel were analyzed on this basis: the blocking effects step in 
immediately as a droplet is formed or (in our case) enters the 
area of the respective segment of the circuit board. The reactant 
supply is interrupted and the current density decreases sharply. 
The performance recovers slowly as an effect of cross-flow from 


channel to channel and cross-diffusion within the gas diffusion 
layer. 


5. Conclusion 


In conclusion, neutron radiography is a unique tool to inves- 
tigate water distribution in unmodified, operating fuel cells. In 
combination with a synchronous measurement of the local elec- 
trochemical performance, not only excess water but also a lack 
of humidification can be detected and much more details caus- 
ing performance fluctuations can be investigated. Variations of 
membrane and gas diffusion layer materials as well as differ- 
ent flowfield designs can be studied and their influence on the 
performance can be investigated which opens a wide field for 
optimization of low temperature fuel cells. 
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